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In a previous paper, based on the black hole perturbation approach, we formulated a new 
analytical method for regularizing the self-force acting on a particle of small mass [i orbiting 
a Schwarzschild black hole of mass M, where fi <C M. In our method, we divide the self- force 
into the S-part and .R-part. All the singular behavior is contained in the .S-part, and hence 
the J?-part is guaranteed to be regular. In this paper, focusing on the case of a scalar-charged 
particle for simplicity, we investigate the precision of both the regularized 5-part and the 
_R-part required for the construction of sufficiently accurate waveforms for almost circular 
inspiral orbits. We calculate the regularized 5*-part for circular orbits to 18th post-Newtonian 
(PN) order and investigate the convergence of the post-Newtonian expansion. We also study 
the convergence of the remaining _R-part in the spherical harmonic expansion. We find that 
a sufficiently accurate Green function can be obtained by keeping the terms up to t = 13. 

§1. Introduction 



We are at the dawn of gravitational wave astronomy. Several ground-based 
interferometric gravitational wave detectors are in various stages of development. 1 ^ 4 ^ 
R&D studies of a space-based gravitational wave observatory project, the Laser 
Interferometer Space Antenna (LISA), 5 - 1 which will observe gravitational waves in 
the mHz-band, are in rapid progress. There is also a proposal for a DECi hertz 
Interferometer Gravitational wave Observatory (DECIGO/BBO). 6) - 7 ) This will be 
a laser interferometer gravitational wave antenna in space sensitive to frequencies 
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f ~ 0.1 Hz, near which all the currently planned detectors are insensitive. 

Among the promising targets for LISA are binary systems consisting of a su- 
permassive black hole (M ~ 10 5_8 Mq) and a compact object of near solar mass 
(/i ~ I-IOMq). Such gravitational wave sources could provide the first high-precision 
tests of general relativity in very strong gravitational regimes. It is expected that 
gravitational waves will be observed from the inspiral stage of these binary systems. 
Because the typical observation time of the inspiral stage is very long, ranging from 
several months to several years, high-order relativistic corrections should be observ- 
able. The standard post-Newtonian approximation 8 ) seems inappropriate for this 
type of binary, i.e. with an extremely large mass ratio, because the extension of the 
PN approximation to higher post-Newtonian orders is not straightforward in this 
case. However, in such a system, there is another natural perturbation parameter, 
/j,/M ~ 10~ 6 . Therefore, the black hole perturbation approach is suited for such a 
binary system. 9 )' 10 ) 

We consider the self-force acting on a particle of small mass \x orbiting a black 
hole of mass M, where [i <C M. The central black hole M fixes the metric of 
the 'background' spacetime. At lowest order in the mass ratio, 0/M)°, the par- 
ticle moves along a geodesic in the background geometry. Already at this lowest 
order, combined with the assumption of adiabatic orbital evolution, this approach 
has proven to be very powerful for evaluating general relativistic corrections to the 
gravitational waveforms. At the next order, however, the spacetime geometry 
is perturbed, and the motion of the particle is affected by its self-force. When we 
consider the point particle limit, the full self- force diverges at the position of the 
particle, and hence it needs to be regularized. The field generated by the parti- 
cle can be expressed in terms of the retarded Green function. The retarded Green 
function can be decomposed into the direct part and the tail part, and the properly 
regularized self-force is given by the tail part (or the -R-part) of the self-field, which 
is obtained by subtracting the direct part (or the 5"-part) from the full field. This 
was first shown by DeWitt and Brehme 11 ) for the scalar and electro-magnetic cases 
decades ago, and rather recently by Mino, Sasaki and Tanaka, 12 ) and Quinn and 
Wald 13 ) for the gravitational case. An equivalent but more elegant decomposition of 
the Green function was proposed by Detweiler and Whiting, 14 ) in which the direct 
part is replaced by the S"-part and the tail part by the i?-part. The 5*-part is defined 
so as to vanish when the two arguments x and x' are timelike, and to satisfy the 
same equation as the retarded Green function. The latter condition implies that the 
i?-part satisfies the source-free homogeneous equation. 

To avoid various technical difficulties associated with the gravitational case, 
many previous papers focused on the regularization problem of the scalar radiation 
reaction force. 15 )~ 25 ) Even in this simplified case of a scalar-charged particle, how- 
ever, only special orbits, such as circular or radial orbits, have been considered, and 
no systematic method for computing the regularized self-force has been given. 

Recently, 27 ^ we formulated a new analytical method for regularizing the self-force 
acting on a particle in the Schwarzschild spacetime. As is commonly the case in the 
self-force regularization problem, our method also uses the spherical harmonic de- 
composition to obtain a regularized expression for the self-force. 15 ^ 22 )' 24 ^ 26 )' 28 - 1-31 ) 
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The novel point of our method is that it provides a new decomposition of the re- 
tarded Green function in the frequency domain, which we call the S and i?-parts. 
This decomposition guarantees that all the singular behavior is contained in the S- 
part. In the black hole perturbation approach, the Green function is conventionally 
calculated in the frequency domain, while the S-part to be subtracted is given in the 
time domain. We presented a systematic method for translating the .S-part into an 
expression in the time domain, and gave an explicit expression for the regularized 
5-part of the self-force, that is, the (S — 5)-part for general orbits. 

Our method for handling divergences in the .S-part relies on the slow motion 
approximation, i.e., the post-Newtonian (PN) expansion. Although the order of the 
expansion is not technically limited, thanks to the systematic calculation method, 33 ) 
the highest PN order attainable in practice is limited by the availability of compu- 
tational resources. In this paper, we revisit the problem of the self-force acting on 
a scalar charge moving in a circular orbit around a Schwarzschild black hole, and 
demonstrate that this practical limitation of our new method is in fact not at all 
severe. We investigate the convergence of the PN expansion and estimate the re- 
quired PN order to obtain sufficiently accurate waveforms. We also clearly elucidate 
the difference between the roles of the conservative part and the dissipative part 
of radiation reaction forces in this formalism. The role played by the conservative 
self-force during the orbital evolution has to this time been demonstrated only in 
toy model scenarios (see, for example, Ref. 32)). We also discuss the convergence of 
the remaining iZ-part. 

The paper is organized as follows. In ^21 we summarize our new regularization 
method in the case of the scalar self-force. In ^31 we calculate the scalar self-force on 
a particle in a circular orbit around a Schwarzschild black hole. First, we explicitly 
obtain an expression for the (S — S) part of the scalar self-force to 4 PN order. It is 
possible to extend the calculation to an arbitrarily high order systematically. In fact, 
though not presented explicitly in this paper, the actual calculation is done to 18PN 
order. Then, using this result, we study the convergence of the (S — 5)-force in the 
PN expansion. Subsequently, combining this result with the calculation of the R- 
part with sufficient accuracy, so that it does not spoil the 18PN order accuracy of the 
(S — (S')-part, the regularized scalar self-force is evaluated and is compared with the 
result obtained by Detweiler, Messaritaki and Whiting. 25 ) In i n order to obtain 
a rough, qualitative estimate of the PN order to which we need to proceed in the 
realistic gravitational case, we pretend that the scalar charge can be replaced by the 
gravitational mass and consider the phase error in the gravitational waveform due 
to the truncation of a series in the PN expansion in our regularization calculation. 
The final section, £JSJ contains conclusion and discussion on the implications of our 
result. 
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§2. Analytic regularization scheme 

We consider a point particle of a scalar charge q moving in a Schwarzschild 
background characterized by mass M, 

ds > = _ U _ IE) dt 2 + (l _ ™L V 1 dr * + r 2 + gin 2 0#2) (2 . 1} 



The equation for the scalar field induced by the particle is given by 



V Q V Q V(z) = -g / rf / 4 *( r )) . (2-2) 



Here, r and z(r) are, respectively, the proper time and trajectory of the particle, 
and ^ is the determinant of the metric. The solution for the (full) scalar field can be 
obtained using the retarded Green function as 



^\ x ) =q dr G tull {x, z{t)) . (2-3) 

J — oo 

The retarded Green function satisfies the Klein-Gordon equation 

V Q V a G full (x, x') = J W ( X - X> ) ; (2 . 4 ) 

with retarded boundary conditions. Furthermore, due to the spherical symmetry and 
the static nature of the background spacetime, the Green function can be decomposed 
in Fourier-spherical harmonics as 

G iull (x,x>) = J ^e— (*-*> Y,d n ±(r,r')Y em (6, fiYSJjf rf) . (2-5) 

£m 

Here the functions Y£ m (0, </>) are the ordinary spherical harmonics. Then, Eq. ()2-4|) 
reduces to an ordinary differential equation for the radial part of the retarded Green 
function, 



2M\ d 2 2(r -M) d ( rJ 2 l{t + 1) 



r / dr 2 r 2 dr \r — 2M r 2 



7 fuU (r _M _ 6 ( T ~ r ') 
A') ' I „2 



r (2-6) 

The radial part of the full Green function is expressed in terms of two independent 
homogeneous solutions of Eq. (|2-6|) as 

^ (r ' r,) = ) WW - r) + ^(r)« " r 1 )) , 

(2-7) 

with Wi mu) being the Wronskian, 



w(«Pin.<Pup) = r I 1 



2M 



A^ p (r)Vr n «-(|:CW)^pW 
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Here (j)\ n is a homogeneous solution with vanishing flux from the past horizon (when 
multiplied by e~ tujt ), and ^>^ p is a homogeneous solution with vanishing flux from the 
past null infinity. They can be obtained using a systematic method developed by 
Mano, Suzuki and Takasugi. 33 ) The parameter v is called the renormalized angular 
momentum, which reduces to t in the limit ujM — > 0. 33 )' 34 - 1 

We can express the homogeneous solutions 4>" n and 4>^ p in terms of another set 
of solutions, (f>% and (p~ u_1 , which are given by a series of Coulomb wave functions 

as 10),33) 



- <t> v c + A, <t>i 



T.C + C" 1 - (2- 



'up 



The properties of the coefficients (3 U and 7^ are discussed in Ref. 33). As shown in 
one of our previous papers, 27 -* using these homogeneous solutions, the Green function 
can be divided into the two parts as 

r') = r') + gf^r, r') , (2-9) 

where where 

gL>,r')= ~* u [€(r)fc U -Hry(r'-r) + <t>-»-\r)r c (ry(r-r')] , 



-1 



?w(r ' ^ (1 - Pv7u)W emu (<l%, fa"' 1 ) 



(2-10) 



The 5-part of the radial Green function, gf mu} , is symmetric, and it satisfies the same 
inhomogeneous equation as the full radial Green function, hence becomes singular 
when the sum over the spherical harmonic indices is taken. Contrastingly, the i?-part, 
9emu> satisfies the source free equation, and hence remains regular. It is the 5-part 
that needs to be regularized. An important fact is that the S-part contains only 
positive integer powers of u when expanded. Therefore the frequency integral can 
be analytically carried out to give the 5-part in the time domain easily. By contrast, 
it is difficult to obtain the i2-part in the time domain analytically for general orbits, 
because it includes terms logarithmic in uj. The contribution of such a term in the 
time domain cannot be expressed in terms of local quantities at the location of the 
particle, but it can be expressed non-locally in terms of the integral along the past 
trajectory. Nevertheless, because the -R-part is regular and the summation over I 
converges, it should be possible to evaluate it numerically without much difficulty. 
The scalar self-force is given by 

F a = qP%d^, (2-11) 

where the projection tensor P a ^ = 5 a ^ + u a u^ is applied to it in order to keep the 
scalar charge constant. The regularized self- force is given by the i?-part of the self- 
field defined by Detweiler and Whiting. 14 ) The i?-part of the self-force in our new 
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decomposition scheme now takes the form 

F* = - F% = (F! - F*) + F* . (2-12) 

Further motivating the above decomposition of the self-force from a physical point 
of view, we note that it has been shown that the (S — 5)-part of the self-force is 
shown to be purely conservative for generic orbits. 27 ) In other words, the radiative 
(or dissipative) part of the force is solely contained in the i?-part. This property is 
an additional merit of our decomposition. 

§3. Convergence test of the PN expansion 

Because our regularization method relies on the PN expansion, it is necessary 
to examine if it converges sufficiently rapidly. First we investigate the convergence 
of a PN series expansion for the (S — «S)-part for circular orbits, as a test of the 
efficiency of our analytic regularization method. To obtain an analytical expression 
for the (S — <S)-part in the time domain for general orbits, it is necessary to expand 
it in powers of oo. This means that we adopt the slow motion approximation, 

ujr = 0(v), loM = 0(v 3 ); v < 1, (3-1) 

where v is the velocity of the particle. A calculation is said to be of n-PN order if it 
is accurate up through 0(v 2n ). 

We next study the convergence of the i?-part. In the circular case, it is easy to 
evaluate even the i?-part analytically in the time domain. For a circular orbit, the 
trajectory z a (r) and the four velocity u a (r) of the particle can be written as 

_ {u'r, r , I, u*r} ; „' = ^Z, n ^ L^ZMZ. . (M) 

Then the frequency integral is readily done by substituting mfi for u>, where Q = 
vft /u 1 . However, unlike the case of the (S — £*)-part, a time-domain expression for 
the .R-part cannot be easily obtained for general orbits. For this reason, we do not 
bother expanding it in u but test the convergence with respect to the summation 
over I. 

In these convergence tests, in order to obtain a rough estimate for the corre- 
sponding tests in the gravitational case, we translate our results in the scalar case 
to the gravitational case by identifying q/VG with the mass /j, of the particle. 

3.1. The (S — S)-part 

The transformation of the 5-part into the time domain makes it possible to 
subtract the divergent 5-part analytically. If we were to perform this subtraction 
numerically, the fraction to be subtracted would become closer and closer to unity 
as £ increases. Apparently, this would imply a stringent requirement of numerical 
accuracy. In this sense, we anticipate a clear advantage of the analytical subtraction. 

As noted above Eq. 1)3-1(1 . it is necessary to expand the (S — 5)-part in powers 
of u, which corresponds to a PN expansion, to obtain its expression in the time 
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domain. Here, however, we would like to emphasize that the PN expansion should 
not be taken as a limitation of our general scheme. Let us explain the reason. The 
decomposition into £ modes is unavoidable, as long as we adopt the mode-sum (mode- 
by-mode) regularization, in which the large I modes are the cause of the divergence. 
Now, once we regularize the 5-part by subtracting off the 5-part, the result is a 
convergent series in i modes, which is essentially a PN series. Thus, we cannot 
completely avoid PN expansion even if we employ a fully numerical approach. 

When the orbit shrinks and approaches the innermost stable circular orbit 
(ISCO), it is believed that the PN expansion is a poor approximation. This be- 
lief comes partly from the expected difficulty in carrying out the PN expansion to 
very high orders, and partly from the possible worsened convergence of the PN ex- 
pansion, particularly for orbits close to the ISCO. We demonstrate that this belief is 
unfounded on both counts. Specifically, with regard to the first point, we calculate 
to 18PN order to show that systematic calculations to very high PN order are indeed 
possible. With regard to the second point, which is the central issue of this paper, 
we find that the convergence is rather rapid, at least for almost circular orbits, and 
the expansion up to 4PN or 6PN order turns out to be sufficient for the (S — »S)-part, 
even for orbits shrinking up to the ISCO. 

The force arising from the (5 — 5)-part is purely conservative for generic geodesic 
orbits. The conservative part of the self- force influences the orbit by pushing a par- 
ticle off the geodesic orbit, and it also appears in source terms for second order 
perturbations. Therefore, it clearly affects the waveforms, but with corrections of 
higher order in ii/M. Specifically, when we investigate the orbital evolution, correc- 
tions due to the conservative self-force are, roughly speaking, given by its ratio to 
the Newtonian force. In contrast, since the radiative part is that responsible for the 
inspiral, errors due to the truncation of the PN expansion at a finite order are to 
be evaluated in comparison with its leading term. Hence, errors in the conservative 
self- force have an additional suppression factor, /i/Af, compared with those in the 
radiative part. This is partly the reason why a rather low PN order, such as 4PN, is 
adequate to obtain the (S — 5)-part with sufficient accuracy. 

The components of the self-force from the (S— 5")-part of the scalar field, F^~ s = 
F^ — F^, have already been obtained for general orbits in Ref. 27). This part of the 
self- force is expressed in terms of local quantities of the particle, i.e., its position and 
velocity. Therefore, all we need is to provide in these general formulae information 
specific to the circular nature of the orbit, Eq. ()3-2j) . Only the r-component is non- 
vanishing for the (S — S')-force, because the t- and (^-components are directly related 
to the rates of change of the energy and angular momentum, and they are purely 
dissipative for circular orbits. It is given explicitly as 



<S-S 




L 



133 



73 




1107284037660637 
400151300487120 



+ -7T 2 ) V 

64 J 




182118981911377689978271 29 vr 2 



) 



V 8 +0{v 9 ) 



(3-3) 



8548630707351386171520 



where V 



\[MI> 



vqQ. Here we have expanded i 7 ^ 



up to 4PN order. However, 
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it is computed up to 18PN order in the convergence test below. 

First, in Fig. ^ we display the convergence of the r-component of the (S — S)- 
force as a function of the PN order for several representative orbital radii Tq. Here, 
an estimator of the convergence of the PN expansion is defined by 



F s- 

± a 



F 



S-S\ 



n-1 



(3-4) 



where F^~ \ n denotes the (S-S)-part of the force truncated at ra-PN order, and the 
denominator F a denotes the exact (fully relativistic) self- force including the R-p&it. 
In practice, because it is impossible to know its exact value, we use the most accurate 
result in our calculation. It is found that the convergence of the PN expansion is 
steady near the ISCO, although it is less rapid there. The convergence improves 
slightly when the Pade approximation is used near the ISCO. Here, in the Pade 
approximation, we have chosen the denominator to be quadratic in V 2 . 

3.2. The R-part and the regularized self-force 

In the case of circular orbits around a Schwarzschild black hole, each component 
of the self-force due to the i?-part is formally given by 



?R 



iq 2 f2 sr , % 



w 



I'm 



F^ — ~T V] dr9itn.mn( r i r 0, 



Cm 



1 pfi 



r=r 



IT 



Y tm (-,0 



7T 



Y im [-,0 



(3-5) 



The i?-force in this case can be completely obtained analytically, because the inte- 
gration with respect to u is done by substituting mQ for u. Also, note that since 
only modes up to a finite value of £ contribute to the self- force for a given PN order, 
the computation of the i?-force is relatively easy. If we need precision up to ra-PN 
order inclusive, it is sufficient for us to calculate the modes up to £ < n + 1. The 
detailed analysis is summarized in Ref. 27). The 4PN results, after summation over 
the £- modes, are 



F 



R 



F 



73 
133 



16151, 
21014 



-V 2 



395567 i 
106808 



-V q 



4 4, . Trt 1196206548879997 \ T . 

-7 + - In (210 + V 

3 ' 3 v ' 400151300487120 / 



+ 



/ 59372120592232147984979 14 



66 



V 1709726141470277234304 



_ i n( y) _ _ in(2) 



14 7 



q z V 
Airrl 



^ + 2 ^- 7 lv^ + %V< i + 0{v" 
3 24 5 



0{v 9 ) 
(3-6 
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Fig. 1. The relative error of post-Newtonian formulas in the r-component of the (S — S')-force at 
ro = 6M, 10M, 20M and 50Af . The horizontal axis is the order of the post-Newtonian expansion. 
The top figure displays the convergence in the Taylor expansion and the bottom figure is that 
in the Pade approximation. 



Here 7 is Euler's constant, 7 = 0.57- • • . The temporal component, which represents 
the energy loss rate, starts at 1.5PN order, as is expected for dipole radiation. 

To see the convergence of both the r-component and t-component of the R 
force, in Fig. El we plot A^(n), defined in a manner analogous to A^~ s (n), for 
several representative orbital radii ro. The i?-part contains the radiative part, whose 
errors are to be compared with its leading term, as mentioned above. Thus we need 
to evaluate the R-p&rt with better accuracy than the (5 — S^-part. In this case, it 
would seem that the convergence of the naive post-Newtonian expansion given above 
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Fig. 2. Plots of errors in the post-Newtonian expansion of both the r-component (upper panel) and 
the ^-component (lower panel) of the 7i-force. 



is too slow for a small r$. Since the -R-part plays no role in the divergences, it is 
not necessary to transform its expression into the time domain analytically. This 
means that we do not have to expand the i?-part of the Green function in powers of 
uj. Here we propose to use just the spherical harmonic expansion (^-expansion) , and 
to compute the contribution from each ^-mode with sufficient accuracy. A criterion 
for truncating various series expansions in g^ nuJ (r, r') is discussed in the Appendix 
in some detail. Plots of the ^-expansion similar to that in Fig. |2] are given in Fig. El 
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Fig. 3. Plots of the errors in the I expansion of both the r-component (upper panel) and the t- 
component (lower panel) of the R- force. Here, A^{1) is defined in the same way as in the case 
of Aa(n), but in the ^-expansion instead of the post-Newtonian expansion. The convergence is 
much faster than that of the naive post-Newtonian expansion. 



3.3. Comparison with Detweiler, Messaritaki and Whiting 25 ) 

The total self-force in the case of a circular orbit around a Schwarzschild black 
hole is obtained by summing the (5 — 5)-part and the i?-part as 



4irr% 



4 7 2 

-—7 H 7T 

3 ' 64 



- 3 H2V)- 2 - 



+ 



/ 604 
Vl5" 



+ 



29vH 
1024 



(3G 



ln(2) 



14 



ln(V) 



14 



-1 



+ 0(v 9 ), (3-7) 
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Table I. The r-component of the self- force at several radii (M = 1 and q 2 = 4ty). 



ro 


6M 


10 M 


20M 


50M 


F?(ro) 


1.676820878 x 1(T 4 


1.378448171 x 10~ b 


4.937905866 x 10~ Y 


6.346791373 x 10~ y 



Ft = F t- (3-8) 

In the r-component of the scalar self-force, there is a significant cancellation between 
the (S — S)-part and the .R-part, and the total force begins at 3PN order. 

Here we compare our result with that obtained by Detweiler, Messaritaki and 
Whiting. 25 ^ They calculated the radial component of the self-force for the case 
ro = WM, M = 1 and q 2 = Air, and evaluated the uncertainty using a Monte Carlo 
simulation. It should be noted that the source term of the field equation in our 
definition is different from theirs by a factor of 4tt. Their result is 

= 1.37844828(2) x 10 -5 . 

On the other hand, with the result obtained using a Pade approximation for the 
(S — S')-force accurate to 18PN order and with the most accurate i?-force in our 
calculation including the terms up to £ = 18, we obtain 

F r R = 1.378448171 x 1(T 5 . 

In Table QJ we list the numerical values for the r-component of the regularized self- 
force for r = 6M (ISCO), WM, 20M and r = 50M. In our computation the 
accuracy is limited by the (S — S')-part. Hence, the accuracy of the full regularized 
force can be read from Fig. ^ 



§4. Errors in gravitational wave cycles 



In this section, we evaluate the PN order necessary to obtain sufficiently accurate 
gravitational waveforms. For this purpose, we consider the number of wave cycles 
N when the particle spirals in from an initial radius rj to a final radius r/. 

For the gravitational wave search with known waveforms, the so-called matched 
filtering method is used. In this method, the correlation between the detector's out- 
put and theoretical templates for the waveforms is taken to search for the maximum 
signal-to- noise ratio (SNR). Because a phase error of order unity in the theoretical 
templates strongly reduces the SNR, an error in the number of cycles, denoted by 
AN, caused by truncating the force at a given order of the expansion is a good indi- 
cator of the significance of the error. Here, because the frequency of the gravitational 
wave is given by twice the orbital period, we estimate twice the number of cycles of 
the orbital rotation by the formula 

i r> n (dEiy 

it J n dE/dt 

"Tk dV (^Wf?' (4 ' 1) 
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where V = y/Mjr i} V f = y/M/r f , and f t R = {4irr%/q 2 ) x F R . 

First, we consider the correction due to the conservative part of the self- force. 
The factor Q(dE / dr®) in Eq. (|4-1|) is determined by the background geodesic motion 
at lowest order. Hence, the self-force contribution is suppressed by its ratio to the 
Newtonian force. This ratio is given in the PN expansion as 



pS-S 

r r _ __H_ \ „ T/ 2i 

(/iM/r 2 ) M 



J>F 2i . ( 4 - 2 ) 



Hence, the correction to N due to the [S — 5)-force at n-PN order is given by 

AN § - S (n) = 4a n dV - (1 ~ ^ R V 2n . (4-3) 
V ) nj^ (i-3V 2 ) 2 f R 

Because the ratio F^~ s / '{pM /r 2 ) contains a factor fi/M, the //-dependence in this 
correction appears only through Vi and Vf. We set Vf equal to the velocity at the 
ISCO, where higher PN order corrections become the largest. Explicitly, Vf is fixed 
to 1/y/E. The value of Vi depends on the masses ix and M and the observation 
period, but here we set Vi equal to in order to estimate the maximum error. This 
makes AN S ~ S (re) completely independent of fi. Table ITU lists the error estimates. 
From this table, it is evident that the expansion up to 4PN order is in most cases 
sufficient for the (S — S')-part of the self-force. However, in the case we expect to 
detect an inspiral signal with a very large SNR, the observational error in the number 
of cycles, AN, may be made small, say, AN < 0.1. In such cases, in order to extract 
as much orbital information as possible, it will be necessary to calculate up to 6PN 
order or higher. 

Next we consider the dissipative part of the self-force, which is completely con- 
tained in the -R-part. In particular, in the case of circular orbits, the t- and ip- 
components of the i?-force are purely dissipative. Because the self-force is responsible 
for the orbital energy loss, the energy loss rate dE/dt which appears in Eq. Q4TJ1 is 
solely determined by F R . Therefore, in contrast to the conservative part, the higher 
PN corrections of the dissipative self- force are not suppressed by a factor of \xjM in 
comparison with its leading order term. Hence, the PN corrections through dE/dt is 
not suppressed by a factor of ijl/M. We therefore define the error indicator AN R {£) 
by 

AN k (£) = N k {t)-N k {t-l), (4-4) 

with 

s M f v f (1 - 6V 2 ) 

N R (£)=4-I dV { — ™! (4-5) 



V JVi (1 - 3V 2 ) 2 f 



Table II. The relative error AN S s (n) for the conservative self-force. 



n 


2 


3 


4 


5 


6 


7 


8 


9 




15.1 


0.660 


0.292 


0.00413 


0.0635 


0.00927 


0.00962 


0.00143 
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where f t is the i?-part of the force summed up to the ^-th harmonics. 

Interestingly, when both the mass of an inspiraling small compact object and 
the observation period are fixed, the error AN R {£) has a maximum as a function of 
the mass of the central black hole. Hence, an upper bound on the required order 
of the ^-expansion can be obtained by evaluating the phase error formula (|4-4|) at 
the maximum. The presence of a maximum can be understood by considering two 
limiting cases: 

• Small M limit: . The largest truncation error comes at the ISCO. Therefore, we 
fix Vf at the ISCO. Then Vi is fixed for a given observation period. Here, we 
adopt 1 year for the observation period. When M is small, Vi is small. Thus 
we can obtain a good upper bound on the error in AN by setting Vi = 0. After 
substituting for Vi, the expression (|4-4|) for AN is manifestly proportional to 
M. Hence, AN decreases as M decreases. 

• Large M limit:. Again, we fix Vf at the ISCO. Then, because M is large, 
the orbit remains close to the ISCO. In the limiting case, we can regard the 
integrand in Eq. (|4-4|) as constant. Then, we have AN ~ NAFt/Ft, and the M 
dependence appears only through N. Because for a given observation period, 
iV decreases as M increases, AN also decreases as M increases. 

Because AN decreases for both small and large M limits, it should have a maxi- 
mum. Therefore, we can identify the "required" order of the ^-expansion for a given 
observation period independently of M. We plot AN R (£) as a function of M in 
Fig. |IJ where we set \i equal to 1 M Q and the observation period to 1 year. The plots 
exhibit a peak near M = 10 5 M Q , and the error decreases as M moves away from 
this value on both sides. From these plots, for an inspiraling compact object of solar 
mass and for a given observation period of O(lyear), the expansion up to £ = 13 is 
found to be sufficient irrespective of the mass of the central supermassive black hole. 
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Fig. 4. Error in the number of cycles caused by truncating the ^-expansion in the i?-force as a 
function of the mass of the central black hole. 
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§5. Conclusion and discussion 

The new regularization method proposed in Ref. 27) is based on the post- 
Newtonian (PN) expansion. To exhibit the effectiveness of our regularization method, 
we have considered the simple case of a scalar charged particle in circular orbits 
around a Schwarzschild black hole, and we have shown here that one can actually 
compute the PN expansion up to sufficiently high orders. 

We have analytically obtained the (S — 5)-part of the scalar self-force up to 
18PN order (though the explicit results are displayed only up to 4PN). Using this 
result, we investigated the convergence of the PN expansion for the (S — S")-part. 
The results are shown in Fig. ^ It appears that our result up to 18PN realizes an 
accuracy of ~ 10~ 7 for tq = WM. At the innermost stable circular orbit (ISCO), the 
convergence slows down, but is found to be steady, and the accuracy of the obtained 
regularized self-force is high, with relative error of O(10 -4 ). 

For the i?-part, we have computed the contribution from each spherical harmonic 
mode up to £ = 18. The expansions with respect to e = 2Muj and U = M/ro are 
truncated at a sufficiently high order. Because the i?-part is unaffected by the 
regularization procedure, we do not have to perform its PN expansion. Interestingly, 
the convergence of the ^-expansion is found to be much faster than the PN expansion, 
although the results truncated at the ^-th harmonics are correct only up to (^+1)-PN 
order. 

Next, we have considered the phase error in the gravitational waveform due to 
truncation of the PN expanded series by interpreting the scalar charge as the mass 
of an orbiting particle. We have found that the 4PN order calculation of the (S — S)- 
part seems to be sufficiently accurate. For the iZ-part, we need a calculation up to 
£ = 13 to make templates for a one-year observation up to the ISCO, assuming that 
the mass of the inspiraling star is 1 Mq . 

Comments are now in order concerning the convergence behavior considered in 
the Appendix. First, we note here that the expression for the Green function given 
in Ref. 33) is convergent for any values of e and z. However, the convergence is 
not guaranteed once we expand the Green function in powers of these parameters. 
In fact, the expansion of the (5 — S")-part of the Green function with respect to e 
has a finite convergence radius. This can be understood as follows. 35 ) In the limit 
e — > 0, we have v — > i. Now, as the value of e increases, v decreases and eventually 
becomes I — 1/2. For this value of u, (j) u and (\>~ v ~ x are no longer independent, 
and the Wronskian W{4> v ', (f)~ u ~ 1 ) vanishes there. Therefore, we have a simple pole 
in 9i mLU ( r > r ')i an d hence the power series expansion fails to converge at this value 
of e. This problem is hidden as far as circular orbits are concerned, because high 
frequency contributions are completely suppressed. However, for generic orbits, we 
have arbitrarily high frequency contributions. For example if we expand the radial 
component of an orbit with respect to the eccentricity e, terms of higher order in 
e have higher frequencies. Hence, uj can be arbitrarily large. These high frequency 
contributions, which cannot be handled by our analytic regularization method, are 
present in general. Nevertheless, we believe that our scheme is useful in a relatively 
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wide region of the orbital parameter space. 

In addition to the problem mentioned above, the spins of the central black hole 
and the inspiraling compact star also shift the orbital frequency at the ISCO. We 
may encounter situations in which relativistic effects become more important than 
we have considered in this paper. Therefore, developing an alternative numerical 
method 36 -* that can play a complimentary role is also important. 

Ultimately, our goal is to derive the gravitational self-force on a point particle 
orbiting a "Kerr" black hole for a generic orbit, and eventually to construct highly 
precise theoretical template waveforms, to be used in the upcoming era of gravita- 
tional wave astronomy. For this purpose, we need to develop the second-order black 
hole perturbation theory. Our present results suggest that a relatively low PN order 
is sufficient for the second-order perturbation. This fact encourages us to develop a 
formalism along the lines of our scheme whose treatment of the singular part exploits 
the PN expansion efficaciously 
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Appendix 

Criteria for Truncating in gf^^ 

Because the infinite sum over the spherical harmonic modes is practically im- 
possible, we need to decide at which value of t to truncate the sum, which should be 
determined by an accuracy requirement. 

For the t-component, if we compute the sum up to I = 18, the accuracy is 
~ 10~ 8 for ro = 6M, as shown in Fig.|3]or Fig. [21 Approximately the same accuracy 
is obtained by the truncation at I = 18 for the r-component as well. However, 
because an accuracy of the r-component is determined by that of the {S — S")-part, 
if we require an accuracy of ~ 10~ 4 for F^~ s , we only need to compute up to £ = 11 
for F k . 

The above conclusion, however, assumes that g^^ for each given t is com- 
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Fig. 5. Plots of the error in the ^-expansion of both the r-component (upper panel) and the t- 
component (lower panel) of _F 1 ' (8,25) , F 2 ' (8,25) and F 3 ' (8,25) (denoted by 1, 2 and 3, respectively) 
at r = 6M. 



puted with sufficiently high accuracy. In this appendix, we explain our method for 
computing gf mijJ and examine its accuracy. Our method is based essentially on the 
post-Minkowski expansion, in which the black hole mass M plays the role of the 
expansion parameter. Specifically, we expand the i?-part of the Green function with 
respect to e = 2Mw = 0(v 3 ) and U = M/r = 0(v 2 ). 

A.l. Analytic expressions for the R-part 

The fi-part of the Green function is given by 



-1 



Mu {€ {r)<t>- v -\r') + 07"" \rW c {r')) 



{l-Mv)W bnu (<%,fc v - x ) _ 



(A-l) 
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Using the variables z = lot and e = 2Mw, the function 0£ and the coefficients (3 and 
7 are given by 



^ n ^r(n + ^ + l + ie) T(2^ + 2) 



-i(z-e) / 2 /„ _ \y 

ej; ^ . -n r , +l + ie) r( 2n + 2v + Z 

n=—oo v ' x ' 

XiFi(n + z^ + 1 + ie; 2n + 2v + 2; 2i(z - e)) , 
r(-t/ + ze) r(2i/ + 2)^_„_! 

r(i/ + i + ic) r(-2i/) ' 
.r(i/ + i + ic) r(-2i/) 



,sin7ru/ + le) 



r{-v + ie) r(2u + 2) sin vr(z^ - ie) ' 
where iiq is the confluent hypergeometric function, and 



(2(z - 6))" 
(A-2) 
(A-3) 

(A-4) 



2-^-^(1 - 2ie) 
r{l + v + ie) 2 r{l + v - ie) 



E 



E 

n=0 



r(n + + 1 + ie)r(n + 2u+l) 



n\r(n + v + 1 — ie) 



(-n)!r(n + 2i^ + 2) 



-i 



We assume u > in the above derivation of the homogeneous solutions. To obtain 
the solutions for uj < 0, we consider their analytic continuation to complex u> through 
the upper half complex w-plane. This is done by replacing to in with |w|e i7r . 
The Wronskian is given by 



2u + 1 

2lo 



E 



(A-5) 



It is important to note that we can easily obtain expressions for a v n (for — N < 
n < N where N is a large integer) to very high accuracy in powers of e, for any value 
of v. In fact, because the convergence condition of the quantities a v n for n — ► ±oo 
determines v, it is the accuracy of v that is important. Therefore, the accuracies 
of the coefficients $, 7 and W£ muJ (4>c , 4>~ u ~ l ) essentially depend only on v. In other 
words, we may regard these as functions of only v when we estimate their accuracies. 



A. 2. Error analysis 

Apart from the dependence on u, which we leave unexpanded in all the quanti- 
ties, we perform the post-Minkowski expansion by replacing z with e/(2U). We then 
assume e = 0(t> 3 ) and U = 0(v 2 ). We express the homogeneous solution cp" as 



Kir). 



(A6) 



We then write a truncated expression for gf mtjJ as 

fl£L Bl,BB) (r, r') = 9 lt un2 \r, r>) + 9 2 ±^\r, r>) + g%£"»Hr, O , (A-7) 
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where 



AM) 
AM) 



v-l 



'r 

(ni,n2) 



("1,^2) 



+ (r <-> r'), 



<I 



(r)*-" _ V) 



(A-8) 



The indices (711,712) indicate that the corresponding quantities are expanded in 
terms of e and v 2 to 0(e 2ni ) and 0(?j 2n2 ), respectively. Note that 711 and 712 cor- 
respond to the post-Minkowski order and post-Newtonian order, respectively. The 
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Fig. 6. The error A" (ri) in v. 
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coefficients A4 are composed of /3, 7, We r , 



e u and U u . As noted above, 



because their accuracy is essentially determined by the accuracy of v, we regard the 
Ai as functions of only v, and we estimate the error caused by using an approximate 
value of v. 

First, we investigate the convergence of v, which is given by a power series in e. 
We define A u (n) by 



A 1 



in = 



}(n) 



i(n-l) 



(A- 



Here, the index (n) indicates that these quantities are expanded in terms of e to 
0(e 2n ). Considering the most relativistic case, we set r = 6M and e = 2M£f2, where 
Q is the orbital angular velocity. We plot A u (n) for various values of £ in Fig. |H1 
We see that the convergence is very fast. Because the errors in Ai are proportional 
to that in u, we may ignore the errors in A® if we take, say, n = 8. 

We now investigate the convergence of the -R-force. Let us denote the part of 
the force that is due to ff^L 1 ( r ' r> ) p^( ni > n2 \ As seen from Fig. |5J the i?-force 
is dominated by p^ ni ' U2 \ This is in agreement with the order-counting carried out 
previously. 27 ^ We define the errors Z\ 2 (rii) and A 2 a (ri2) by 



Al{m) = 



p2,(ni,n 2 ) _ p2,(ni-l,n 2 ) 
"a F a 



pR 
^ a 

j-,2, (711,712) T-,2,(ni ,712—1) 
F a — F a 



pR 



(A-10) 



The t-components of Z\ 2 (?t,i) at n<i = 25 and A 2 a { n 2) at n\ = 8 are plotted in Fig.[7| 
and Fig. |HJ respectively. We see very fast convergence. Although not shown here, 
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we find that the same is true for the other components. These results confirm that 
our calculations are sufficiently accurate. 
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